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The disagreement of the scahng of the correlation length ^ between experiment and the Ginzburg- 
Landau (GL) model for domain chaos was resolved. The Swift-Hohenberg (SH) domain-chaos model 
was integrated numerically to acquire test images to study the effect of a finite image-size on the 
extraction of ^ from the structure factor (SF). The finite image size had a significant effect on the SF 
determined with the Fourier-transform (FT) method. The maximum entropy method (MEM) was 
able to overcome this finite image-size problem and produced fairly accurate SFs for the relatively 
small image sizes provided by experiments. 

Correlation lengths often have been determined from the second moment of the SF of chaotic 
patterns because the functional form of the SF is not known. Integration of several test functions 
provided analytic results indicating that this may not be a reliable method of extracting ^. For both 
a Gaussian and a squared SH form, the correlation length ^ = 1/cr, determined from the variance 
of the SF, has the same dependence on the control parameter e as the length ^ contained explicitly 
in the functional forms. However, for the SH and the Lorentzian forms we find ^ ~ 5^^^. 

Results for ^ determined from new experimental data by fitting the functional forms directly to 
the experimental SF yielded ^ ~ £~" with u ^ 1/A for all four functions in the case of the FT 
method, but v ~ 1/2, in agreement with the GL prediction, in the the case of the MEM. Over 
a wide range of e and wave number fc, the experimental SFs collapsed onto a unique curve when 
appropriately scaled by ^. 

PACS numbers: 47.54.-|-r,47.52.+j,47.32.-y 



I. INTRODUCTION 

Spatially extended nonlinear non-equilibrium systems 
continue to be of great interest because they yield quali- 
tatively different phenomena that do not occur in linear 
systems 0. One of these phenomena is spatio-temporal 
chaos. The examples of spatio-temporal chaos found 
in Rayleigh-Benard convection (RBC) lend themselves 
to particularly detailed experimental study under excep- 
tionally well controlled external conditions 12] • RBC oc- 
curs in a thin horizontal layer of fluid with thickness d 
heated from below when the temperature difference AT 
exceeds a critical value ATc 3]. 

Particularly noteworthy is the state that occurs in 
RBC when the fluid layer of density p and shear viscosity 
•q is rotated about a vertical axis with angular frequency 
Lo. When the dimensionless frequency O = iocPp/rj ex- 
ceeds a critical value, then the pattern immediately above 
onset consists of disordered domains of convection rolls 
known as domain chaos. This is illustrated in Fig.^ and 
^p. Within each domain the roll orientation is more or 
less uniform; but different domains have different orien- 
tations HmSimilllSini- The domains are unsta- 
ble and undergo a persistent but irregular dynamics. In 
this case, the time-averaged mean-square velocities and 
temperature deviations from the conduction state grow 
continuously from zero as ATc is exceeded, i.e. the bi- 




FIG. 1: AU images are for Q. = 17.7 and e = 0.05. Upper 
left: shadowgraph image of size 60d x 60d from F — 61.5. 
A movie is available 12]. Upper right: shadowgraph image 
of size 44d x 44d from F = 36. Lower left: the average SF 
S {kx, ky) computed using the FT method and a Kaiser-Bessel 
window with a — 2.5 and averaged over 4096 images like the 
one in the upper left. Lower right: the average SF from same 
data as shown in lower left, except S{kx,ky) was computed 
with the MEM. 
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furcation is supercritical. 

Chaos immediately above a supercritical bifurcation 
offers a unique opportunity for theoretical study because 
weakly-nonlinear theories are expected to be applicable. 
These theories, in the form of Ginzburg-Landau (GL) or 
Swift-Hohenberg (SH) equations, by virtue of their gen- 
eral structure and from their numerical solutions, pre- 
dict that the inverse half-width at half-height ^ of the 
structure factor (SF, the power spectrum of the pattern) 
should vary as e"" with v = 1/2 s.s e = AT/ATc - 1 
vanishes El El El El 113, 111 El , and that it is the 
only length scale in the problem. Thus it was particu- 
larly disappointing that measurements for domain chaos 
disagreed with this expectation Q. Determinations of 
a correlation length ^ based on the variance of the SF 
found f oc e^'^'ff with Veff — 1/4 rather than 1/2. 

Several possible explanations of the apparent disagree- 
ment with theory were explored by various authors. Hu 
et al. proposed that defects and fronts injected by the 
side wall into the bulk may play an important role |llj |. 
Laveder et al. demonstrated that additive white noise 
intended to mimic the effect of wall defects decreases v 
to the value of v^f / measured experimentally provided 
that the noise level is sufhciently large '2^. Recent un- 
published experiments in our group using a sample with 
a radially ramped spacing j23, H^l render the side-wall- 
injected defect idea an unlikely explanation. On the basis 
of numerical simulation using a SH equation Cross et al. 
proposed that the finite size of the experimental convec- 
tion sample decreases the effective value of v . 

Experimental determinations of a characteristic length 
scale ^ usually are based on numerical estimates of the 
variance cr^ of the SF, with 1 = l/a H ll^. This is so 
because the precise analytic form of S{k) is not known in 
the nonlinear regime above onset. In Sect.^we examine 
the relationship between ^ and the length <^ that appears 
explicitly in various functional forms that might be used 
as approximations to the SF. We find that ^ ~ S}/"^ when 
^ is the length appearing in the SF of the linear GL or 
SH equation. Thus, in those cases ^ does not have the 
same e-dependence as ^. However, we find that ^ ~ ^ for 
the square of the linear SH SF and for a Gaussian form 
for the SF. The results suggest that the success of the 
moment method of data analysis depends upon the rate 
at which the SF decreases at large k. 

In Sect, mil we discuss two methods of determining 
the SF from patterns. The standard method, using the 
Fourier power spectrum, is sensitive to the limited size 
of the experimental images l25!|. To overcome this limi- 
tation we also employed the maximum entropy method 
(MEM), which is discussed in detail in that section. The 
benefits of the MEM are exhibited through analysis of im- 
ages from a SH simulation detailed in Sect. II Vl We find 
that the MEM is quite powerful in its ability to overcome 
the finite image-size problem. 

In Sect.0we present new experimental results for pat- 
terns in the domain-chaos state. We compare results 
from the Fourier analysis to the results from the MEM. 



We determined ^ by fitting several possible functional 
forms for the SF to the data and found that the results 
for ^ are not very sensitive to the form of the fitting func- 
tion. In the case of the Fourier analysis, the fits did not 
change the result i/ ~ 1/4 that had been found before 
by the moment method. The reason why the moment 
method also gave this result can be found in the behav- 
ior of S{k) at large fc, where S{k) drops off fast enough 
for ^ to be essentially proportional to ^. In the case of the 
MEM, we found that ~ 1/2 as expected. This indicates 
that the findings from Fourier analysis are dominated by 
the finite image-size effect. 

We also examined the maximum height B of the SF 
and found B ^ with /3 ~ 3/4 in the case of Fourier 
analysis. The two results v and /? ~ 3/4 conspire 

to retain the expected dependence of the total power on 
6"+'^ with V + P = 1 even though both v and [3 do not 
have the expected value 1/2. Once again, we find that 
the MEM overcomes the finite data length yielding (i ~ 
1/2. We also find that the result for v and (3 are not 
strongly dependent on f2 for either the Fourier analysis 
or the MEM. Thus we conclude that, in light of careful 
analysis using the MEM, the experimental domain-chaos 
state does possess a length scale in agreement with the 
prediction from the GL model. 

Finally, in Sect. IVII we examine the extent to which 
the SF can be represented by a unique scaling function 
over a range of e and k. We obtain excellent collapse 
of the data, both in the case of Fourier analysis and the 
MEM, when the results for ^ obtained from a fit to the 
SF at each e are used to scale the SF. 



II. MOMENTS OF ANALYTIC FUNCTIONS 

In order to better understand the consequences of using 
numerical moments to compute the correlation length, we 
derived analytic expressions for the zeroth, first, and sec- 
ond moments of several proposed forms of the SF. This 
yielded results for the correlation length based on mo- 
ments that could be compared with the correlation length 
^ that occurs directly in the functional form chosen for 
the SF. We found that the ^ ^ ^™ where m depended on 
the particular form of the structure factor. This suggests 
that using moments to measure the correlation length 
does not necessarily yield with the e dependence im- 
plied by the GL model. 

We investigated four particular forms of the SF. The 
exact form for domain chaos is not known. For our pur- 
poses one useful form is the SF of the linear SH equation 



e{k^~klf + Akl 

It is an excellent approximation to the SF of the lin- 
earized full equations of motion (the Boussinesq equa- 
tions) in the presence of additive noise for RBC below 
onset |2^ . Another useful form is the squared SH SF 
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S{k) = 



(2) 



used by some authors for fits to numerical results in order 
to estimate a half width at half height 5k of the distribu- 
tion [27l,l28l|. In the limit of large ^, the correlation length 
^ in Eq. n approaches l/5k, but in Eq. |21 approaches 
l/(5fc for large 6- Whe n we directly compare ^ 
and as in Fig. ^1 below, we divide by Vv2 — 1 so 
that it approaches 1/Sk. 

We also considered a Gaussian form 



in the limit of large ^. In the case of the squared SH 
form, 



P = 2nklB 



which reduces to 
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S{k)^Bexp - {ln2) {k - kof S,^ 
and a Lorentzian form 



Sik) 



B 



eik-koY 



(3) 



(4) 



The latter is the SF of the one-dimensional linearized GL 
equation. The correlation length ^ of both the Gaussian 
and the Lorentzian is exactly equal to l/5k regardless of 
the size of ^. For all four of these functions we note that 
the position of the peak of S (k) is at fco and S (fco) = B. 

The total power is P = 27r S (fc) kdk. In the case 
of the SH form, 



in the limit of large ^g. 

The first moment is fc = (27r/F) S (fe) k'^dk. 
In the case of the SH form 



k = ko 
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which reduces to 



P = 



(6) 



in the limit of large ^. 

In the case of the squared SH form 



k = 




(11) 



which reduces to 



fco 



(12) 



fc2 



^ 25, fco 


f-tan-1 (-^)' 
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in the limit of large . 

The second moment does not converge for the SH form 
but it does for the squared SH form. In the case of 
squared SH, cr^ = fc^ 



fc , where 



and Eq. ^2 gives fc . For large ^, 
1 



es 



1 - 
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so that ^ ^s- ond moment remained finite. In that case cr^ 

In order to compute a similar expression for the 2'kP~^ S {k)k^dk — k so that 
SH form, we introduced a cutoff kc so that the sec- 
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FIG. 2: 5 computed from integration of the SH SF and the 
squared SH SF. The analytic results were evaluated with ^ = 
1.6e"^/^ fc = TT, and kc = 14. Solid line: C = 1.6e"^/^ 
(chosen to match the solid line in Fig. lib! . Dashed line: ^ from 
SH computed with Eq. 1151 using kc = 14 and computing k 
numerically so that the integration is done only over the range 
< fc < 14 for consistency. Dashed-dotted line: Approximate 
5 from SH computed with Eg. 1161 Dotted line: computed 
from Eq. 1111 and Eq. 1131 Thin solid line that passes almost 
directly through dotted line: Approximate from squared 
SH computed from Eq. ^] 



and in the limit of large ^ 




Equation ^| indicates that in the case of the SH form 
C ~ -s/C regardless of the cutoff kc- Although we do not 
present the details here, we also found that ^ ~ \/J for 
the Lorentzian form but ^ ^ ^ for the Gaussian form. 

Figure 121 compares these results for ^ computed by nu- 
merically integrating the SH SF and the squared SH SF. 
The lowest order expansion for the squared SH SF in 
Eq. 1141 is an excellent approximation to the exact re- 
sult as it lands almost directly on top of it. However, 
although the exponent is —1/2, it does give a that is 
smaller than the input ^ (shown as the solid line in Fig. [21) 



because of the extra constant \/ \/2 — 1 discussed earlier. 
The approximate result for the SH SF (from Eq. Ilt)|) is 
not very close to the exact result (from Eq. I15|l because 
the exact result for k (Eq.jSJ was used to compute Ea.lTBl 
while the dashed curve shown in Fig. |21 was computed 
with a finite cutoff for k to match the cutoff used in the 
integral for fc^ for the sake of consistency. Nevertheless, 
the approximate result is parallel to the exact result with 
cutoff indicating that regardless of these details, the ex- 
ponent from integrating the SH SF is —1/4. 

According to the GL model, at small e one expects 
^ ~ e-" with 1/ = 1/2 [H |l7|. One also expects the 
total power to vanish at onset in proportion to e. Thus, 
if Eq. ^or Eq. [Ogives the shape of 5* (fc) correctly, then 
according to Eq.j^lor Eq.jSl P ~ e'^+^with f + p = 1 and 
B - e'^ with P = 1/2. 



III. ESTIMATION OF 5* (fc) FROM PATTERN 
IMAGES 

We found that the accuracy of S{k), estimated from 
the Fourier power spectrum, suffered greatly due to the 
finite size of the images of patterns that are available from 
experiment. In applying the Fourier-transform (FT) 
method to compute S (fc) from experimental data, we 
divided the images by a reference image taken below on- 
set, applied a Kaiser-Bessel window ^9] to the divided 
images, and computed the magnitude squared of the FT. 
The azimuthal average of the squared magnitude of the 
FT yielded the SF S'(fc). 

We experimented with three windowing functions: a 
square window, a Welch window, and a Kaiser-Bessel 
window |29j 

(i) ^ I lo (TTa^l-it/rf^ /lo (™) \t\ < r 

\ \t\>T. 

(17) 

In Eq. El t is the distance along one of the axes from 
the center, r is the half-width of the image, and Iq is 
the modified Bessel function of the first kind and order 
zero. The parameter a controls the rate at which the 
window drops off as the edge of the image is approached. 



5 



Since the data was two dimensional, the window was the 
product of a one-dimensional window function in each 
direction. Since the windowing function attenuates the 
signal around the edges of the image, it reduces the total 
power. To compensate for this we divided S (fc) by a 
constant so that the total power in the final SF agreed 
with the total power of the raw image. 

The specific windowing function did not greatly affect 
the result for ^. We ultimately settled on using a Kaiser- 
Bessel window with a = 2.5 for the results given in this 
paper. Our major results, namely the exponents v and /? 
and the scaling collapse of the SF, were independent of 
the windowing function. 

We used a fast Fourier-transform (FFT) algorithm js^l 
capable of transforming images of arbitrary size to com- 
pute the SF so that no interpolation or zero padding was 
required to re-size the image to an integer power of two. 
We found that interpolation distorts the large-fc behavior 
of the SF because it smooths out random noise that is 
present in the original signal. The use of zero padding 
circumvents the smoothing of the noise and thus avoids 
this distortion at large k, but instead inserts ringing at 
small k. Due to the easy availability of FFT algorithms 
that can work on any size image and the speed of modern 
computers, there is no reason to sacrifice the behavior of 
the SF at either large or small k. A minor tradeoff is that 
previous work utilized an aspect-ratio correction to the 
ima ge t hat removed a slight anisotropy in the frame grab- 
ber |3ll| . Because this correction requires interpolation of 
the raw image we avoided using it. We expected only a 
very slight error by doing so because the raw images are 
nearly square. However, in the F = 36 sample there was 
a radial distortion, due to optical aberration, that was 
strong enough to warrant its removal in order to avoid a 
systematic error in the length scale. We did not inves- 
tigate the large k behavior in that sample, so we expect 
no significant problem from the distortion correction. 

Even after dividing the experimental images by an op- 
tical background, the resulting SF still contained the non- 
deterministic part of the background spectrum present in 
the images taken below onset. Since fluctuations |33.l33| 
are too feeble to be detected for the parameters of the 
present experiment, we attributed this background sig- 
nal primarily to electronic noise and subtracted a back- 
ground, determined below onset, from the SF above on- 
set. As a result, for the case of experimental data, the 
analysis that follows is applied to the background sub- 
tracted SF 5S {k) EE S{k)- Sb (fc) where Sb (k) is the 
SF averaged over many images below the onset. We also 
applied the FT to simulations of the SH domain-chaos 
model. In that case there was no need to divide by a ref- 
erence image or to subtract electronic noise, so we used 
S (fc) computed directly from the simulation images. 

We attempted to overcome the finite image-size prob- 
lem of the FT method by using the maximum entropy 
method (MEM) to estimate S (fc). Although this method 
is commonly used for ID data |34j . computing the spec- 
trum of 2D data with this technique is still somewhat 



of an open problem 1351 . We implemented the algo- 
rithm detailed in Ref. [36j which uses an iterative method 
to arrive at the power-spectrum estimate. Figure 4 in 
Ref. provides a detailed flowchart of the MEM algo- 
rithm, which we followed precisely. The MEM provides 
the power spectrum as an expansion of the form 



S {k'Xl ky) 



1 



(18) 



where A (ni, 77.2) are the coefficients of the expansion and 
F [■ ■ ■] is the discrete Fourier transform. Because the 
MEM provides an expansion containing sines and cosines 
in the denominator, in contrast with the FT where sines 
and cosines are in the numerator, the resulting power 
spectrum may more accurately represent sharp peaks be- 
cause it may contain poles l34| . It was also shown previ- 
ously [33 that the MEM is much less sensitive to short 
data lengths than standard Fourier analysis. 

The algorithm in Ref. 36j depends on an accurate es- 
timate of the auto-correlation of the data. A central sec- 
tion is cut from the auto-correlation data and used in 
the iterative process. For all the analysis in the present 
work we used a region size of 45 x 45 data points, which 
corresponded to about 12d x 12d for the F = 61.5 sample 
and also for the SH simulation. The iterative process pro- 
duces auto-correlation data that is continued beyond this 
region according to the spectral-entropy maximization- 
criteria that defines the MEM. The number of coefficients 
A (ni, 7^2) determines the size of the continued region. In 
principle it is possible to use a relatively small continued- 
region size, in order to reduce the required amount of 
computation, and then embed the resulting A (711,712) 
into a larger region, setting higher order coefhcients to 
zero in order to produce a finely meshed power spec- 
trum. In practice we found that this approach, while 
advocated in Ref. js^, did not always produce positive- 
definite power spectra and thus was not stable for our 
purposes. Instead we ran the entire algorithm on a large 
720 X 720 mesh for A (711, 712). Although this was heav- 
ily computation intensive, it did provide a reliable power 
spectrum estimate after a sufficient number of iterations. 
On our modest fleet of 2 GHz PowerMac G5s, we could 
run roughly 200 iterations per minute per CPU. Typically 
we needed 50,000 to 1 million iterations for convergence, 
depending on the e step, making the processing of all the 
data in the present work a fairly massive undertaking. 

The number of iterations depended on the auto- 
correlation error eo, as defined by Eq. 25 in Ref. l36|, 
which was reduced to a specified level. In the case of the 
SH simulation data discussed in Sect. II Vl we found that 
eo = 10~^ was sufficiently small to give convergence of the 
SF. Table m shows the results for 1/ and /3 in the conver- 
gence test. By eo = 10^^, v and /3 have nearly reached a 
constant. In the case of the experimental data discussed 
in Sect. EI we found that we only needed eo = 10~^ for 
satisfactory convergence. This is likely due to the bet- 
ter statistics for the auto-correlation function in the case 
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TABLE I: Convergence test for SH simulation data with F* 
150. 



10" 
10~ 
10" 
10" 



V 

0.858 
0.619 
0.555 
0.533 



^.211 
0.455 
0.496 
0.506 



of the experiment due to the fact that we averaged the 
auto-correlation function over 4096 images per e step in 
the experiment, but only 256 images per e step in the 
simulation. 

We applied the MEM to the pattern images in a sim- 
ilar way that the FT was used. We averaged the auto- 
correlation over many images in order to reduce statis- 
tical fluctuations. Then the MEM was applied to this 
averaged auto-correlation to yield a single estimate of 
S {kx,ky) at each temperature step. This is in contrast 
to the Fourier analysis where we computed an S {k^, ky) 
for each image, directly from the Fourier power spec- 
trum, and then averaged them together to get an aver- 
aged spectrum. Since the Fourier transform is linear, this 
technique is algebraically equivalent to making a single 



Fourier transform of an averaged auto-correlation func- 
tion, similar to what was done for the MEM case. Due 
to the iterative nature of the MEM algorithm we used, it 
was prohibitively slow to compute an S {k^, ky) for each 
image individually. 

We used the SF to determine the correlation length. 
However, we modified the analysis scheme compared to 
earlier work ^ |2J| , in light of the results presented in 
Sect, m to avoid using numerical moments. We de- 
termined ^ by fitting each of the functions in Eqs. ^ 
21 multiplied by the shadowgraph transfer function from 
Refs. ^,^1, to the data near the peak of the SF aver- 
aged over all images at a particular e. 



IV. FINITE IMAGE SIZE EFFECT IN 
SIMULATION OF SH MODEL FOR DOMAIN 
CHAOS 

The finite size of the pattern images has a significant 
impact on the accuracy of a measurement of ^ from S (k). 
To study this effect, we ran simulations using the Swift- 
Hohenberg model for domain chaos |3 . We utilized the 
algorithm in Ref. ^3 ^.nd periodic boundary conditions 
to solve the equation 



J 



I 



(19) 



for ip, a field that can be used to model the temperature 
of the convection sample at the midplane. FigureOshows 
an example of ip. The time step for numerical integration 
was 0.1. The initial condition for ijj was a random grid of 
straight-roll patches. At each e, 10000 warmup time steps 
were performed followed by 256 snapshots of tp recorded 
at an interval of 1250 time steps. The pixel spacing was 
chosen to reflect a non-unity sample thickness, unlike the 
choice in Ref. (17i |. in order for the wave number and ^ 
to be nearer the values in the experiment. The choice of 
this constant has no effect on the value of (3, or /i, thus 
it does not affect our main conclusions. 

The control parameter in the simulation is related to 
the experimental control parameter by e = ^ — 

2.60e m, where kc is the critical wave number and is 
the curvature of the neutral curve [40j . The numerical 
value 2.60 corresponds to 17 = 17.5. Figure^lshows both 

and as a function of as computed numerically 

from the neutral curve. We note that the exact value of 
effect on the results for ^, B, f3, /, or 
/i, and only serves to adjust the range of e. 

The parameters gi, g2, and 53 control the stability 
balloon and can be chosen to model a specific f2 and 
Kiippers-Lortz angle Okl- In the present work we used 



gi — 1, g2 — —2.4534, and 53 = 0.522 which are the same 
parameters used in Fig. lb of Ref. ^3 ^.nd which corre- 
spond to 9k L — 51° and roughly fi ~ 17.5 for the Prandtl 
number in the present experiment. Note that Eg. 1191 as- 
sumes infinite Prandtl number, so this rough value of 
comes from comparing il/f2c, where flc is Prandtl num- 
ber dependent. 

As a check on the accuracy of our implementation of 
the solver algorithm we attempted to reproduce Fig. 7a 
of Ref. 17]. Our simulation agreed within 2% for the 
data points at e — 0.1, 0.2, and 0.3 in that figure, but 
differed by about 27% for the data point at e = 0.03 for 
an unknown reason. Nevertheless, we are confident in the 
correctness of our simulation and in the results to follow. 

In order to determine the effect of the finite image size 
on the measured value of ^ we ran simulations of image 
size 1024 x 1024 corresponding to an image aspect ra- 
tio r* = 300, where F* refers to the horizontal width 
of the image. This is in contrast to the aspect ratio F 
which refers to the physical aspect ratio of the convec- 
tion sample. From these large simulations we cut various 
sized center sections and computed the SF using both 
the FT and the MEM following the procedure described 
in Sect. IIIl( We fit the SH function and the squared SH 
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FIG. 3; A solution of Eq.mfor e = 0.12. This is a 512 x 512 
(r* = 150) cutout from the center of a 1024 x 1024 image 
with r* = 300. A movie is available |4H . 
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FIG. 4; Top figure: curvature of the neutral curve. Bottom 
figure: 4/fc2^g. 
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FIG. 5: Structure factor (SF) from simulations using Eg. 1191 
with e = 0.12, computed with the Fourier-transform method. 
Both data sets are from the same simulation except with a 
different size center-cutout used to compute the SF. Triangles: 
F* = 300. Circles: F* = 75. Solid lines: Fits of Eq.[T]to the 
data over the range fc ± 3/^. 



function to the SF. 

The choice of F* strongly affected the shape of the 
SF when computed with the FT. Figure [S] compares the 
SF of the same data for two different values of F* . The 
same data analyzed with the MEM, as shown in Fig. El 
is affected slightly by F* , but it is not nearly as sensitive 
as the FT. It is important to note that the sharpest SF 
peak is not necessarily the best. The F* = 300 SF from 
the FT is sharper than either of the MEM peaks, but it is 
not at all consistent with its shorter data length relative, 
F* = 75. In contrast, the MEM SFs land nearly on top of 
each other when comparing F* = 300 and F* — 75. The 
ability to accurately represent the SF peak for relatively 
small image sizes is crucial because system sizes as large 
as F* = 300 are experimentally inaccessible. 

From fitting the SF for both the FT and MEM, we 
computed ^ for many e and F* as shown in Fig. [7| Once 
again we found that the MEM provides significantly more 
consistent results for vastly different F* . In the case of the 
FT results, it is critical to note that increasing F* does 
not simply increase ^ due to the larger data length. It also 
increases the slope of ^ vs. e on logarithmic scales, thus 
affecting the measured v. Since we wish to measure v 
accurately enough to compare with the model prediction 
V = 1/2, this is a discouraging outcome. 

We ascertained the effect of image size on the accuracy 
of V by computing v s& & function of F* as shown in 
Fig-IHl The two lowest F* points in that figure correspond 
to the available image size in the experimental samples 
discussed in the present work. The v values from the 
MEM and FT method both approached nearly the same 
asymptotic value at large F* . Fortunately, the result from 
the MEM has nearly reached this asymptote at even the 
smallest F* indicating that it may be reliably used for 
measuring v even in the case of small images. Not only 
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FIG. 6: Structure factor (SF) from simulations using Eq. 1191 
with e = 0.12, computed with the maximum entropy method 
(MEM). Both data sets are from the same simulation except 
with a different size center-cutout used to compute the SF. 
Triangles: F* = 300. Circles: F* = 75. Solid lines: Fits of 
Eq.0to the data over the range k ± 3/^. 
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FIG. 8: Results for v from fits like those shown in Fig. |7| 
Solid symbols: from the MEM. Open symbols: from the FT 
method. Circles: ^ was determined by fitting the SH SF Eq.0 
to the data. Squares: ^ was determined by fitting a squared 
SH SF Eq.|21to the data. Dashed line: the v — 1/2 prediction 
from model equations. 
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FIG. 7: Measurements of ^ from fits of Eq.^to many SFs like 
those shown in Figs. |S] and [H] Open symbols: from the FT 
method. Solid symbols: from the MEM. Triangles: F* = 300. 
Circles: F* = 75. Solid lines: Fits of the power law ^ oc e~'^ 
to the data over the range shown. 



is the FT measurement of v much too low at the smallest 
r* , but it also contains rather large fluctuations at those 
values. At the largest F* it is almost as good as the MEM, 
but that does not help for the analysis of experiments 
which are limited roughly to F* < 100 ^3 ■ 

We observed a similar finite image-size effect for the 
scaling of B with e. As discussed in Sect. ^ we expect 
that B ^ £^1"^ provided that ^ obeys the scaling in the 



amplitude model, i.e. ^ ~ e 



-1/2 



Our SF fits that yielded 



^ also gave B. Figure El shows some examples of B for 
the same data as shown in Fig [3 

We extracted a scaling exponent /? by fitting the equa- 
tion i? cx to the data. As indicated by Fig. IIUI the 



FIG. 9: Resuhs for B from fits of Eq. [H to many SFs like 
those shown in Figs. |5] and |n] Open symbols: from the FT 
method. Solid symbols: from the MEM. Triangles: F* = 300. 
Circles: F* = 75. Solid lines: Fits of the power law B oc e'^ 
to the data over the same range as shown in Fig. Q 



result is similar to the finite image-size effect on v. At 
the smallest F* values, the MEM results exhibit a slight 
F* dependence, but they quickly reach a nearly constant 
value as F* increases. The FT method suffers from a 
strong F* dependence. However, it is quite satisfying 
that for F* near the experimental values, we obtained 
ly ~ 0.25 and P ~ 0.75 from the FT method, in agree- 
ment with what we (and others previously) observed in 
the experiment. 

FigurelTTlshows the results for i^ + fS for both the MEM 
and the FT method. Both are close to the expected re- 
sult 1/ + f3 = 1 , and approach it even more closely as F* 
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FIG. 10: Results for /? from fits like those shown in Fig. |U] 
Solid symbols: from the MEM. Open symbols: from the FT 
method. Circles: ^ was determined by fitting the SH SF Eq.0 
to the data. Squares: ^ was determined by fitting the squared 
SH SF Eq.Hto the data. Dashed line: the /9 = 1/2 prediction 
from model equations. 
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FIG. 12: Determinations of the characteristic frequency scale 
/ of the domain chaos from the auto-correlations of time-angle 
plots of the SFs computed with the FT method. Triangles: 
F* = 300. Squares: F* = 150. Circles: F* = 75. Crosses: 
F* = 58.6. Pluses: F* = 42.2. Solid line: Fit of the power law 
/ (X e** to the data over the same range as shown in Fig. Q 
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FIG. 11: Values of 1/ + P from combining the results shown in 
Figs. 151 andllUI Solid symbols: from the MEM. Open symbols: 
from the FT method. Circles: ^ and B determined by fitting 
the SH SF Eq. 0to the SFs. Squares: ^ and B determined 
by fitting the squared SH SF Eq.Elto the SFs. Dashed line: 



increases. 

Unlike for ^ and B, the finite image size does not in- 
terfere with the determination of the characteristic fre- 
quency scale / of the domain chaos (the domain pre- 
cession frequency) when the FT method is used to com- 
pute the SF. We computed / using the method explained 
in Ref. We extracted a time scale from the auto- 

correlation of the angle-time plot of the radially averaged 
SF. This yielded the results for / shown in Fig.^l All of 
the / data in that figure land nearly on top of each other 
regardless of F*. Since the measurement of / depends 



FIG. 13: Measurements of fi from fits, like those shown in 
Fig-El to data for / computed from SFs with the FT method. 



more on the position of the peak in the two-dimensional 
SF than on its width, the FT method is able provide 
enough resolution to accurately measure / at smaller F*. 

We fit a power law / cx to / because the GL model 
predicts such a power law with /i = 1 . Figure 1131 shows 
that this exponent is independent of F*. However, the 
result fj, ~ 0.9 is somewhat lower that the predicted value. 
We do not know the reason for this. However, the value 
we found for fj, is much larger than the experimental value 
/i ~ 0.6 m. Recently it was found 0| that this difference 
is caused by the centrifugal force which influences the 
experiment but is omitted in Eq. 1191 
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FIG. 14: Fits of Bans. [Till to the experimental SF computed 
with the FT for Q = 17.7 and e = 0.05. Sohd line: fit of SH 
SF Eq.IH Dashed line: fit of squared SH SF Eq.H Dashed- 
dotted line: fit of Gaussian SF Eq. |3 Dotted line: fit of 
Lorentzian SF Eq. 0] All fitting functions are multiplied by 
the shadowgraph transfer function. 



V. NEW EXPERIMENTAL MEASUREMENTS 

We acquired new data using an apparatus described 
previously |38| |. There were two samples which both 
used compressed SFg gas. The primary sample was at a 
pressure of 20.00 bars and mean temperature of 38.00°C 
where the Prandtl number rj/ (pn) {n is the thermal dif- 
fusivity) was 0.87. The aspect ratio T = D/ {2d) {D is 
the sample diameter) was 61.5. A secondary sample with 
F = 36, which is discussed briefly in this work, was pres- 
surized at 12.34 bars with a mean temperature of 38.00°C 
and the Prandtl number was 0.82. 

We determined ^ by fitting each of the functions in 
Eas. lll4l multiplied by the shadowgraph transfer function 
from Refs. |38Ll3^ . to the data near the peak of the SF. 
Typically we used a range fco ± 3/^ where the initial guess 
of ^ and fco for the fitting range was estimated by Eq. 
using the zeroth numerical moment, peak value, and peak 
position. The fit was not very sensitive to the range of 
k provided that the range included a sufficient number 
of points. Figure ITU shows that for the FT method, all 
of the functions provide a good fit near the half height, 
where ^ is determined, and near the peak where the fits 
give the value of i?. Figure^lshows similarly good fitting 
results for the SF estimated with the MEM. 

FigurelTHIshows the effect of the fitting function on the 
dependence of ^ on e for = 17.7. Although each of the 
fitting functions gave slightly different values for ^, all of 
them gave nearly the same value for v, with v ~ 0.25 from 
the FT method and ly ~ 0.5 from the MEM. This suggests 
that the function fitting is a robust method for measuring 
^ because knowledge of the exact functional form of the 
SF is not required to yield a consistent measurement for 
5. In spite of this consistency for v, we concluded that 
some fitting functions are better than others as will be 
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FIG. 15: Fits of Egns. lTlll to the experimental SF computed 
with the MEM for Q = 17.7 and e = 0.05. Top figure: linear- 
linear plot. Bottom figure: log-linear plot. Solid line: fit 
of SH SF Eq. El Dashed line: fit of squared SH SF Eq. H 
Dashed-dotted line: fit of Gaussian SF Eq.|3] Dotted line: fit 
of Lorentzian SF Eq. |1] All fitting functions are multiplied 
by the shadowgraph transfer function. 



seen below. 

The results from the MEM are in much better agree- 
ment, than the results from the FT method, with the 
prediction 1/ — 1/2 from the GL model. This difference 
is even more dramatic when considering a smaller-sized 
cutout region from the data. Figure [T71 shows ^ for two 
different F, where the analysis has been applied to a cen- 
tral region of identical size in either case. Although in 
the case of F = 61.5 the data shown in Fig.^lis from the 
same raw data as shown in Fig. 1161 the FT method gives 
a somewhat smaller value of v because of the reduced 
image size of the smaller cutout region. The MEM also 
suffers some minor decrease in the value of v, as this 
cutout size corresponds to the smallest F* data point in 
Fig.|Sl Although only barely noticeable in the FT result, 
the MEM clearly shows that <^ is larger in the case of 
F = 36 as compared to F = 61.5. A likely explanation 
is the increased influence of the centrifugal force in the 
larger sample. The presence of the centrifugal force has 
been shown to decrease the domain size j42j . 

The diminishing slope at very small e shown in Figs. ll6| 
and 1171 may be indicative of a physical finite-size effect 
f23i | due to the lateral boundary of the experimental sam- 
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FIG. 16: Measurements of (, by fitting Eans. lTlll to tiie exper- 
imental SF for Q — 17.7. Open symbols: SF was computed 
using a Kaiser-Bessel window with a = 2.5. Solid symbols: 
SF was computed with the MEM. Dashed line: e"^''". Solid 
line: e'^^^. Circles: ^ from fitting of SH SF. Squares: ^ 
from fitting of squared SH SF. Diamonds: ^ from fitting of 
Gaussian SF. Triangles: ^ from fitting of Lorentzian SF. 
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FIG. 17: Comparison of ^ for = 17.7 and both F = 36 and 
r = 61.5, measured with SFs computed from 44ci x lid square 
cutouts using both the MEM and the FT method. Squares: 
r = 36. Circles: F = 61.5. Open symbols: FT. Solid symbols: 
MEM. For F = 36 the fitting range was 0.018 < e < 0.11, 
MEM: u = 0.57, FT: u = 0.17. For F = 61.5 the fitting range 
was 0.017 < e < 0.10, MEM: ly = 0.45, FT: u = 0.15. 



FIG. 18: e"^/^ scaled against ^ and F for SI = 17.7 and both 
r = 36 and F = 61.5, measured with SFs computed from 
44ci X 44ci square cutouts using the MEM. Squares: F = 36. 
Circles: F = 61.5. The data is the same as in Fig. 1171 and is 
shown over the range of the fits in that figure. 
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FIG. 19: Total power for Q = 17.7 from SFs computed 
with both the FT method and the MEM and given by var- 
ious methods. Solid line: . Open symbols: from the FT 
method. Solid symbols: from the MEM. Diamonds: power 
from zeroth numerical moment of the SF data integrated over 
ko — 3/^ < A: < fco -|- 3/^. Circles: power given by Eq.|5]using 
the parameters acquired from fitting the SH SF to the data. 
Squares: power given by Eq.|7|using the parameters acquired 
from fitting the squared SH SF to the data. 



pies. However, over the range of fitting shown in Fig. 1171 
no such effect is observed. As a comparison with the re- 
sult in Fig. 3 of Ref. Fig. [H shows a plot of £-1/^ 
scaled with ^ against e~^^'^ scaled with F. The data is 
roughly constant over the range of e indicating agree- 
ment with the GL prediction for i^. Unfortunately the 
data from two different F do not collapse onto a single 
curve as in Ref. . This may be due to the effect of the 
centrifugal force, which was neglected in that work, but 
is unavoidable in the present experiment. 



Although, in the case of the FT method, the scaling 
exponent ly did not agree with the GL model, the scaling 
of the total power agreed for both the FT method and 
the MEM. Figure Uni shows the total power for n = 17.7 
measured by several methods. The zeroth numerical mo- 
ment of the data did not perfectly agree with the total 
power computed from Eq. or Eq. [3 but it was close. 
In the case of those equations, the fit parameters B, fco, 
and ^ (or £_s for the squared SH) were used to calcu- 
late P. Since Eqs. [3 and [7\ represent the total power 
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FIG. 20: B determined by fitting the FT experimental SF to 
the SH function. Solid lines: power law fits to the data to 
measure /3. Pluses: Q. = 15. Circles: Q. = 16.25. Squares: 
Q. = 17.7. Triangles: ^ = 19.5. Diamonds: Q. = 21.7. 
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FIG. 21: ^ determined by fitting the FT experimental SF to 
the SH function. Solid lines: power-law fits to the data to 
determine v. Pluses: = 15. Circles: Q, = 16.25. Squares: 
Q. = 17.7. Triangles: = 19.5. Diamonds: Q. = 21.7. 



over the range < fc < c» and the zeroth numerical 
moment was necessarily comptited over a finite range 
fco — 3/^ < k < ko + 3/^, the slight discrepancy in total 
power is not surprising. A numerical integration of Eqs.^ 
and 121 over a finite k range yields a total power closer to 
the numerical moment result, however, it is absent from 
Fig. El because it is too close to the other data points 
to be easily distinguishable. In spite of this minor dis- 
crepancy, all three measurements of the total power are 
proportional to e as expected. 

Since, in the case of the FT method, the dependence 
of f on e differed from the GL prediction of e-^/^, it 
is important to also investigate the dependence of B on 
£. As shown in Fig. E| the combination of B and ^ 
according to Eqs. and |S1 yielded a total power that 
depended on e in the predicted way for both the FT 
method and the MEM. 

Figure EUl shows the result for B measured by fitting 
the SH SF to the experimental SF at all H. values. The 
slight variation of the values of /3 (slopes of the lines in 
the figure) did not depend systematically on Q and most 
likely it is due to experimental error. The power-law fits 
shown in the figure were over the same range as was used 
to determine v in Fig. One sees that B alone did 
not depend on s^^"^ as expected; instead B ^ with, 
averaging over the results from all 51, /3 ~ 0.73 while the 
r2-averaged v ~ 0.25. In other words, in the case of the 
FT method, B and ^ conspired to produce the expected 
P e^+f^ with 1. 

Figures 1211 and 1231 show the analogous results from the 
MEM. The MEM yielded much more steeply sloped f 
vs. e curves as indicated by Fig. (221 The resulting 
was in much better agreement with the GL model pre- 
diction. Likewise, (3 was much closer to 1/2 than for the 
FT method. Averaged over v ~ 0.46 and /3 ~ 0.63 
from the MEM. FigureEHsummarizes the behavior of the 
exponents for both the MEM and the FT method as a 
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FIG. 22: B determined by fitting the MEM experimental SF 
to the SH function. Solid lines: power-law fits to the data to 
determine (3. Pluses: Q — 15. Circles: Q — 16.25. Squares: 
Q = 17.7. Triangles: n = 19.5. Diamonds: Q = 21.7. 

function of fi. The MEM clearly gave the closest result in 
agreement with the prediction v = 1/2 and also /3 = 1/2. 
However, the FT method yielded the best agreement with 
the prediction v + P = 1, with the MEM not much fur- 
ther off. Considering the results of Sect. lIVI the MEM is 
more reliable. Thus we conclude that the experimental 
length scale is in agreement with the prediction from the 
GL model. 



VI. SCALING OF 5S (k) FOR EXPERIMENTAL 
DATA AND ANALYTIC FUNCTIONS 

The scaling analogy to critical phenomena extends be- 
yond the power-law dependence of f and i? on e. In 
order to more deeply probe this scaling, we re-scaled 
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FIG. 23: (, determined by fitting the MEM experimental SF 
to the SH function. Solid lines; power-law fits to the data to 
determine v. Pluses: SI = 15. Circles: f2 = 16.25. Squares: 
Q = 17.7. Triangles: n = 19.5. Diamonds: = 21.7. 



FIG. 25: The normalized SFs for Q, = 17.7 computed with 
the MEM. Open circles: e = 0.017. Solid circles: e = 0.025. 
Open triangles: e — 0.033. Solid triangles: e — 0.05. Open 
squares: e = 0.08. Solid squares: e — 0.10. 




evaluating the integral normalization over the range < 
k < oo, all parameters cancel, leaving only a function of 
X. This yields 



S{x) = 



1 



27r2(l + 2;2) 



in the case of the SH form, 



S{x) 



TT^{1+X^y 

in the case of the squared SH form, and 



(20) 



(21) 



FIG. 24: Exponents from data shown in Figs. I20I23I Open 
symbols: from the FT method. Solid symbols: from the 
MEM. Circles: ly. Squares: /3. Diamonds: u + p. Solid 
line: u = (3 = 1/2 prediction from GL model. Dashed line: 
iy + P=l. 



the SFs at different e in an attempt to collapsed them 
all onto a single curve. We followed the procedure of 
Ref. fis^. First we normalized the structure factor so 
that 27r SS{k)kdk = 1. In applying this normaliza- 
tion to the experimental data, we evaluated the integral 
over the range fco — 3/^ < k < fco + 3/^. This included 
most of the total power present in the experimental data. 
Figure l^ni shows the MEM results for several e values and 
demonstrates that normalization alone is not sufficient to 
collapse the data onto a unique curve. It is also necessary 
to re-scale the SF on both the abscissa and ordinate axes 
so that cc = (fc - fco) C and SS (x) = {k/S) SS {x). 

Applying this variable transformation to the analytic 
forms in Eqs. I1I4I provides some insight into the effect 
of rescaling. In the critical-point limit, of ^ — > cx) and 



Six) 



exp {—x"^) 
27r3/2 



(22) 



in the case of the Gaussian form. We refrain from list- 
ing the rescaled Lorentzian because the zeroth moment 
diverges, making it impossible to normalize without in- 
troducing a cutoff. 

Figure 1^ shows the re-scaled experimental SFs from 
the MEM on a linear (top figure) and on a logarithmic 
(bottom figure) scale. Although not shown, we found 
comparable results for both the normalized and rescaled 
SFs computed with the FT method. The figure also in- 
cludes the rescaled SH and squared SH SFs given by 
Eos. EHl and [TTI We omit the Gaussian form in Eq. 122 
in order to maximize the clarity of the figure. It is the 
poorest empirical representation out of Eqs. I20I22I and 
the least physically justifiable. Note that the curves of 
Eqs. 1201 and [21 are not fits to the experimental data; 
there are no free parameters for fitting. None of them 
provide a particularly good scaled representation of the 
data near the peak. The experimental results fall some- 
where between the SH (dotted line) and the squared SH 
(dashed line) form. 
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FIG. 26: Scaled SFs for ^ = 17.7 computed with the MEM. 
Top figure: hnear-Unear plot. Bottom figure: log-linear plot. 
The symbols are as in Fig. 1251 Dotted line: rescaled SH 
SF Eq. EOl Dashed line: rescaled squared SH SF Eq. |2l] 
Solid line: numerically rescaled SH SF using finite limits for 
its normalization. Dashed-dotted line: numerically rescaled 
squared SH SF using finite limits for its normalization. 



The disagreement between the scaled data and the 
scaled functions is primarily due to the normalization 
of the experimental data, which is over a finite range un- 
like the normalization in Eas . 1201 and 12 II which is over all 
k. Also shown are collapsed forms of Eqs. ^ and |21 com- 
puted numerically using the normalization over the same 
finite range as in the experimental SFs. The numerically- 
computed forms were evaluated at finite e, instead of 
taking the limit of ^ — > cx), as was done in the analytic 
cases of Eqs. I2UI22I As a result, ^ and feg remained as 
constants in the collapsed forms. We used the values 
of ^ and fco from the e — 0.10 data point. This choice 
barely affected the numerically collapsed curves because 
fco^ was relatively large for all values of e shown. These 
curves provide insight into the agreement of the shape 
of Eqs. n and [5] with the experimental data. The nu- 
merically computed SH SF agreed quite well with the 
collapsed experimental data, while the squared SH SF 
did not. 



VII. SUMMARY AND CONCLUSION 

In this paper we carefully examined two methods, the 
FT method and the MEM [s^l, for determining critical 
parameters describing the spatio-temporal chaos state in 
rotating Rayleigh-Benard convection known as domain 
chaos. We found that a correlation length ^, equal to the 
inverse half-width at half-height of the structure factor 
S{k), can be determined reliably by fitting one of several 
trial functions to the experimental data. We examined a 
Lorentzian, Swift-Hohenberg, squared Swift-Hohenberg, 
and Gaussian form. We prefer this fitting method over 
computing numerical moments from the data (as had 
been done in the past) because we were able to show 
analytically that the length f = l/cr (cr^ is the variance 
of the data) will be proportional to ^ only when S{k) 
decreases sufficiently fast at large k. 

For a technique of estimating the SF, we found that 
the MEM was much superior to the classic FT method. 
There is a severe dependence on image size in the case of 
the FT method that is quite apparent from our analysis 
of various sized images from the simulation using the SH 
domain-chaos model. Since the experimental images are 
relatively small, the MEM is required in order to accu- 
rately determine the SF from the patterns. 

We analyzed new experimental shadowgraph images of 
domain chaos for a sample of aspect ratio F = 61.5 over 
the range 15 J; < 22. Using the above four functional 
forms as fitting functions, we obtained results for ^ that 
depended only slightly on the function used, but that all 
had the same dependence on the distance e from the on- 
set of convection. In the case of the FT method, we found 
^ ~ e~'^ with V ~ 0.25, in agreement with previous mea- 
surements for a sample with F = 40 but in disagreement 
with expectations based on the weakly-nonlinear ampli- 
tude model. Fortunately the MEM was able to overcome 
the image-size problem of the FT method and yielded 
V ~ 0.46, roughly in agreement with the theoretical mod- 
els. We also determined the maximum height B of S{k) 
from fits of the functions to the data and found that 
B ^ £^ with f} ~ 0.73 for the FT method, in disagree- 
ment with the value j3 — 1/2 suggested by the amplitude 
model. However, we note that the FT method yielded a 
total power proportional to e'^^^ with v + (3 ^ 0.98, that 
agrees with the theoretical expectation that v + (3 ~ 1. 
The MEM yielded (3 ~ 0.63, somewhat closer to the pre- 
diction of /3 = 1/2 than the FT method. 

We also showed that it is possible to present the struc- 
ture factor in a scaled form that largely collapses the 
data onto a unique curve. This was the case for both 
the FT method and the MEM, even though the scaling 
exponents of the FT method results did not agree with 
the prediction from theory. 

In summary, our new experimental data and analy- 
sis yielded results that differ from the earlier result that 
u ~ 0.25 in disagreement with theory. We have shown 
that this difference is due to the limitations of the FT 
analysis method. The MEM is capable of extracting seal- 
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ing exponents from the data that are in agreement with 
the theoretical prediction. In addition to the exponent 
V, we also examined a scaling parameter /3 that describes 
the height of the structure factor, and found it to also 
agree with predictions when the MEM results were used. 
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